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GRB 030406 - an extremely hard burst 
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Abstract. Using the IBIS Compton mode, the INTEGRAL satellite is able to detect and localize bright and hard GRBs, 
which happen outside of the nominal INTEGRAL field of view. We have developed a method of analyzing such INTEGRAL 
data to obtain the burst location and spectra. We present the results for the case of GRB 030406. The burst is localized with the 
Compton events, and the location is consistent with the previous Interplanetary Network position. A spectral analysis is possible 
by detailed modeling of the detector response for such a far off-axis source with an offset of 36.9 °. The average spectrum of 
the burst is extremely hard: the photon index above 400 keV is -1.7, with no evidence of a break up to 1.1 MeV at a 90% 
confidence level. 
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1. Introduction 

The INTEGRAL satellite detects gamma-ray bursts (GRBs) 
in two different ways: for a small number of events 
that f all in the field o f view (FoV) of the imager 
IBIS (Uberti nTeTail 12003) and of the spectrometer SPI 



Vedrenne et al. 2003), INTEGRAL provides accurate posi- 
tions (~2 arcmi n) for rapid ground- an d space-based follow- 
up observations (Mereg hetti et al.l2003l) . A significantly larger 
number of GRBs occurs outside of the FoV of the two 
instruments. These bursts can be monitored by the SPI 
anti-coincidenc e system: SPI- ACS dvon Kienlin et alJ |2003; 
iRau et al.ll2005l) . Up to now, the localization of the bursts that 
happened outside of t he field of vie w was only possible with 
the aid of the 3 rd IPN faurlevl 19971) . 

In this paper we show that for some of these bursts it is 
also possible to perform a more detailed localization analysis 
using the Compton mode of IBIS, provided that the burst is 
sufficiently strong and spectrally hard. 
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We present the capabilities of the Compton mode using 
an example of GRB 030406, a burst that was detected out- 
side of the INTEGRAL field of view. GRB 030406 has been 
detected by SPI-ACS on-board INTEG RAL as wel l as by 
Ulysses, Konus, and Mars Odyssey (Hurlev et al. 2003). It was 
a long burst lasting ~ 65 seconds. Ulysses reported the fluence 
of 1.3 X 10~ 5 ergcrrT 2 (25-100 keV). The observations by the 
three satellites yielded an estimate of the position centered on 
RA(2000) = 19'' l m 43.0 s , DEC(2000) = -68° 4' 39.4", with 
the 3cr error-box of 77 square arcminutes. 

In section |2 we describe the Compton mode of IBIS 
INTEGRAL and demonstrate its localization capabilities. In 
section 12.21 we present the spectral analysis of GRB 030406 
while section|3]contains the discussion. 

2. Observations and data analysis 

The IBIS telescope is an imaging instrum ent on-board the 
INTEGRAL satellite with a coded mask ( Goldwurm etalJ 
20031). There are two dete ction layers i n the detector plan e 
ISGRI JLebrun et alJl2003l) and PICsIT JLabanti et al.ll2003h . 
ISGRI is an array (128 x 128) of pixels made of semiconduc- 
tive CdTe, sensitive to photons between 15keV and ~1 MeV. 
ISGRI works in photon-by-photon mode. PICsIT having the 
same detection area as ISGRI, is an array of 64 x 64 Csl scintil- 
lators, sensitive between ~ 1 70 ke V and 1 5 MeV located 94 mm 
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below ISGRI. PICsIT (usually) makes light curves of the whole 
detector by making histograms of all events (so called Spectral 
Timing Mode). ISGRI and PICsIT can act as a Compton tele- 
scope, registering photons that are scattered in one and ab- 
sorbed in the other detector. The coincidence time window is 
a parameter programmable on-board and it was set to ~ 4 mi- 
croseconds at the time of GRB 030406. The Compton mode is 
sensitive between 200 ke V and ~5 Me V. 

There is a 3 meter long collimator on top of the detector 
unit. The walls of the collimator are made of lead and act as a 
shield to photons with energies up to » 200 ke V. For the geo- 
metrical reasons the optical depth of the shield is smaller for 
photons arriving at large angles. Thus hard photons from off- 
axis sources can pass through the shield and reach the detector. 
In particular off axis GRBs may be detectable in IBIS. 

Thanks to the Compton mode data we were able to make 
successful, independent burst localization. Including informa- 
tion from the ISGRI data we have performed burst spectral 
analysis. 

2.1. Compton imaging and localization 

We have analyzed the Compton mode count rate during the 
GRB 030406. The burst profile in the Compton mode fol- 
lows clearly the profiles seen by SPI-ACS, ISGRI and PICsIT 
(Fig. Q. The gamma ray burst was detected in the Compton 
mode at the level of 30cr reaching ~400 ctss -1 in the sec- 
ond peak (average background amount to ~100 ctss -1 ). The 
Compton mode telemetry did not suffer the data gaps in con- 
trast to ISGRI. 

The Compton mode provides us with the following infor- 
mation: 

- energy deposits in ISGRI (Ei) and PICsIT (£ P ), 

- position of the detection in ISGRI (xy, yi) and PICsIT (xp, yp), 

- timing of the event. 

The two positions provide the information about the di- 
rection of the scattered photon. The two energies provide the 
Compton scatter angle 6c, which cosine, in case of the forward 
scatter, is given by: 



cos 8c — 1 - 



m e <r 



E\ E\ + Ep 

where m e c 2 is the electron rest energy. Given the angle 0c and 
the direction after scattering we obtain a ring in the sky which 
contains the possible directions of the primary photon. In case 
of a point source all the rings cross at the location of the source. 
For a given observation a collection of such rings provides a 
Compton map of the sky S . In real observations the Compton 
map S consists of the source signal and the background B. 

In the case of GRB 030406 we selected the Compton events 
from the begin of the peak up to the end of the tail, which 
consist of a 9.75 s time interval, see Fig.^ For these events a 
Compton map S y of the sky was calculated using the method 
described above with the angular resolution of one degree. In 
order to determine the background we used the pre-burst data 
divided into 200 intervals lasting 9.75 s each. For each interval 
we calculated the all-sky Compton map and we determined the 
mean background B,j and its variance Vtj in each pixel. 
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Fig. 1. Lightcurve of GRB 030406 seen by the INTEGRAL 
sub-systems: SPI-ACS, IBIS/ISGRI, IBIS Compton and 
IBIS/PICsIT The gray areas: precursor, peak and tail corre- 
spond to the time regions chosen for spectral analysis. Time = 
corresponds to the SPI-ACS trigger time at UTC 2003-04-06, 
22:42:03.23. 

The map <x /; was obtained by dividing the difference be- 
tween the Compton map S u and the background B (/ by the map 
of the square root of variance of the noise in each pixel JVij. 
We present the map of significance <x (i in Fig. [2] The largest de- 
viations from the noise are ~20, and the position of the pixels 
with the strongest deviation is consistent with the IPN location. 

In order to obtain an estimate of the accuracy of such posi- 
tioning we used two approaches: the maximum likelihood and 
the Monte-Carlo. In the first approach we assumed that the 
probability density that the true position lies in a given pixel 
i, j is proportional to the likelihood function: 



IBIS 
PICSIT SPTI 




£ u oc exp (cr, 7 ) . 



After normalizing such probability density we found the region 
containing 68% of the probability corresponding to the ~ 3.5° 
radius. 

In the se cond approach we used the IBIS mass model 
(Laurent et a l. 2003) to extract the signal from simulated GRB 
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Fig. 2. All-sky Compton map taken during the second peak of the GRB 030406. The inset shows the dashed area around the 
GRB 030406 centered on the IPN position. 
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Fig. 3. The normalized cumulative distribution of the offset of 
best fit localizations of simulated bursts obtained in the Monte 
Carlo simulation: 68% of the events are localized within ~ 4° 
of the actual position. 



events inserted into stretches of data with experimental back- 
ground. We then performed the identical data analysis proce- 
dure and found the localization of such simulated bursts. We 
have analyzed 127 simulated bursts, and found the cumula- 
tive distribution of the distance between the actual burst po- 
sition and the estimate obtained in our analysis. This distribu- 
tion is presented in Fig. [3] The Icr region (containing 68% of 
the cases), corresponds to the accuracy of ~4°, and the 2cr is 
constrained to 6.5°. The two methods lead to consistent results. 



2.2. Spectral analysis 

In addition to the Compton mode GRB 030406 has been also 
detected by ISGRI and PICsIT There are only two energy 
channels in the PICsIT data stream therefore we have decided 
to model together the Compton and ISGRI data, neglecting the 
PICsIT. 

We have divided the data into three time intervals marked in 
gray in Figure[Oas: precursor, peak and tail. The division of the 
main peak of the burst into the peak and tail was a consequence 
of the ISGRI data loss during the burst due to telemetry gap. 

We first ran Monte Carlo simulations using the mass model 
of the INTEGRAL spacecraft to generate the detector response 
matrices (DRM) for a source located at the position of the burst. 
The DRM matrices were obtained for the ISGRI and Compton 
sub-systems. For all of them we have taken into account the 
influence of the IBIS VETO system, the energy thresholds in 
each ISGRI pixel and PICsIT crystal, the ISGRI noisy pix- 
els, the failed ISGRI pixels and PICsIT crystals as well as the 
ISGRI and PICsIT energy resolution. 

The Compton data were binned into six energy channels 
covering the range from 200 to 2500 ke V. The ISGRI data spans 
from 24 to 450 keV which partially overlaps the Compton en- 
ergy range and were binned into seven channels. The analysis 
of the pre-burst data with the correction for the varying dead 
time effect and exposures provided an estimate of the back- 
ground. This background was subtracted from the total count 
spectrum in each channel. 

We used the standard XSPEC 1 1 .2 JArnaudll996h to fit the 
data. For each time region we tried to fit several models: sin- 
gle power law, broken power law, and a non physical model 
of single black-body. The single power law model never fitted 
the combined ISGRI and Compton data. The broken power law 
fits the data very well in all analyzed regions: precursor, peak 
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Fig. 4. GRB 030406 spectrum fits to the combined ISGRI and Compton mode data. Upper panels show count spectra, while lower 
corresponding vF v spectra. Broken power law model was fitted in each case: to the precursor data (left), to the peak part data 
(center) and to the tail (right). See text for details 



and tail. In additional fit we have allowed for different normal- 
ization of ISGRI and Compton components for peak region, 
taking to account two facts related with dead time influence: 

- during precursor and tail the count rate grows with factor 
about 1 .5 for ISGRI and Compton mode, while during peak, 
these factors reach 2.5 for ISGRI and 4 for Compton mode, 

- the width of the burst peak, «ls, is much shorter than the 
INTEGRAL dead time effect calculation bin, equal 8s. 
Except meaningful decreasing of the reduced x 2 we have not 
find significant evidence of changing fitted parameters. Results 
of all fits are presented in Tabled The broken power law fits 
and the residuals are presented in the top panels of Figure|4] 

We note that the precursor and the tail region can also be 
fitted with black body spectra with temperatures w lOOkeV. 
This is due to the fact that at high energies the count rates fall 
steeply. The lower panels of Figure 0] show the deconvolved 
vF v spectra in each time interval. The precursor and tail spectra 
have peak power around 300-500 keV, while the vF v spectrum 
in the peak seems to rise all the way into high energies. 

The high energy index of yS = -1.7 + 0.3 for peak may 
imply the absence of a peak in the source power spectrum vF v 
in the range covered by our data, i.e. up to 2.5 MeV. In order 



to estimate the lower limit on such a break compatible with the 
data we conducted the following analysis. We fit the Compton 
data with a broken power law with fixed high energy index: 



f(E) = N 



*-* "^ *-* break 



break 



(i) 



t^ •> tL break 

In this initial fit the break energy was constant 
£break=3 000keV. This value was chosen high enough not 
to influence the power law part of the fit. For this fit we 
obtained x 2 - 0.55 with 4 degrees of freedom. We used the 
standard model parameter estimation dLampton et alJ 1 1976) 
to calculate the values of x 1 corresponding to 90% and 99% 
confidence regions, which are respectively Ax 2 — 8.3 and 
Ax 2 = 13.8. Then keeping all the parameters but Ebreak 
constant we looked for the lowest value of £break for which 
X 2 was equal to the above values. We present the results in 
Figure |5] We obtained £break >1 HOkeV at 90% confidence 
level and £break >880keV at 99% confidence level. 

3. Discussion 

We have shown that the IBIS Compton mode is able to de- 
tect and localize GRBs outside of the field of view of the 
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Table 1. Results of the spectrum fits for three time ranges of the GRB 030406. All the errors are at the lcr level. 
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Fig. 5. Confidence limits on ^cutoff : tne solid line represents the 
X 2 as a function of iicutoff using the model of equation Q. We 
present the x 2 values corresponding to the 90% and 99% con- 
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INTEGRAL telescopes. We have found that the position of 
GRB 030406 obtained using the Compton mode technique is 
consistent with the IPN localization. The position accuracy has 
been estimated using the maximum likelihood and the Monte 
Carlo method. The results of both methods are consistent and 
the lcr errors are ~ 4°. 

Given the position of the burst we computed the DRM for 
a source outside of the field of view of the satellite at this par- 
ticular location. Using the DRM we analyzed the the Compton 
data, along with the ISGRI detection and fitted the spectrum 
with a broken power law in the range from 50 ke V to 3 Me V. 

The peak spectrum is very hard; in the vF v the low energy 
index below 400 keV rises with the index a +3.5 and above 
this energy it is still positive » +0.3. We note that this spec- 
trum is an average over the time interval of 2.8s around the 
maximum, thus the spectrum of the maximum might have even 
been harder since the burst spectra usually evolve from hard in 
the peaks to soft in their tails. We do not see any evidence that 
the vF v spectrum peaks below 1 . 1 MeV. Thus, there is a hint 
of existence of bursts with th e peak of vF v spectrum above the 
distribution shown bv lMallozzi et alJ Jl995) which truncates at 
~ 1 MeV for the bright BATSE bursts. 

Such hard spe ctra have already been seen by BATSE - 
Pre ece et al.l (12000) present the distribution of the low and high 
energy spectral indexes for a sample of 156 bright bursts. Both, 
the low and the high energy spectral indexes of GRB 030406 



lie in the upper tail s of the distrib utions of the spectral indexes 
presented bv lPreece et alJ 1120001) . It is therefore clear that the 
spectrum of this particular burst is in clear contradiction with 
the synchrotron model of GRBs (Katz 1994; Tavani 1995), 

which predicts that there is a strict upper limit on the low en- 

7 J T ^J. 

ergy spectral index of -j (Preece et al. 1998). The low energy 

spectral slope is consistent within the error bars with the jitter 
synchrotron model of GRBs (Medvedev 2000). We note that a 
detailed study of time resolved spectra (Ghirland a et alJl2003l) 
showed that the low energy spectral slopes are large ~ 0.5 - 1 
at the rising parts of several bursts. 

The value of E p of GRB 030406 - the peak of the vF v spec- 
trum lies above 1 . 1 MeV. It should be noted that thi s is on the 
high end of the Amati and Ghirlanda relations l l Amati et alJ 
20021 iGhirlanda et aljEooH) . If this relation holds and we as- 
sume that the peak energy is larger than the 90% lower limit of 

1 lOOkeV, then the isotropic radiated energy for GRB 030406 is 
larger than (6 - 20) ■ 10 53 erg, where the uncertainty stems from 
the inaccuracy of the Amati relation. Hard bursts are therefore 
very energetic and may come from a population residing at high 
redshifts. 

In summary, we confirm that the INTEGRAL in the 
Compton mode can detect hard GRBs. The localization accu- 
racy is a few degrees. The GRB spectra can be studied in the 
range from ~200keV to ~3 MeV depending on the location of 
the burst in the instrument coordinates. Further study of such 
hard bursts detectable in the INTEGRAL Compton mode may 
yield potential candidates of high redshift burst and also test the 
limits of validity if the Amati relation. Our rough estimations 
of the Compton mode capabilities give 2-3 detections of such 
bursts per year. They are consistent with BATSE results and 
seem to be confirmed by the other detections (Marcinkowski et 
al. in preparation). Finally, we note that analysis of Compton 
data in principle allows to determine the polarization. A study 
of of the polarization sensitivity of the IBIS Compton mode is 
currently under way (Laurent et al. in preparation). 
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